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Abstract 

The solid-on solid (SOS) model in two dimensions (d = 2) is now solved under the 
constraint of constant energy and then under the new constraint of constant total area. 
From the combinatorial factors g(E; L, M), the new ensemble is constructed with its free 
energy F(A tot ,T) of a membrane of constant (onedimensional) area A tot . The entropy 
per column Y = (1/L) log g(E; L, M) of rectangular L x M strips reduces to a common 
curve in reduced variables. Definitions of the "area" and of the interfacial tension, are 
compared or discussed. Analytical calculations are supported with numerical ones and 
vice versa. Overall the constraint reduces the ratio of A tot to the projected area L, as 
compared with canonical calculation, strongly at high temperatures. 
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I. Introduction 

The solid-on-solid (SOS) model originated very long time ago in connection with 
crystal surfaces and crystal growth; its first extensive description is found in the clas- 
sical paper of Burton, Cabrera, and Frank [1]. It was later realized that the model of 
irregularities of a crystal face could be fashioned into a theory of morphology and also 
of fluctuations of surfaces. Then it was discovered that the SOS model of equilibrium 
surface could be solved exactly in two dimensions, i.e. for a one-dimensional interface, 
by the method of transfer matrix [2-7]. There are reviews available [8- 10]. The solution 
could be applied to free energy calculations and to the theory of capillary waves [7, 11, 12] 
of interfaces. 

Here we apply the SOS model to the difficult problem of membrane fluctuations. 
All interfaces and all membranes as well, undergo shape fluctuations. These necessarily 
result from thermal motion and are intertwined with all thermal fluctuations. However, 
there are important distinctions between interfaces and membranes. An interface is an 
open system, also open with respect to particle exchange. A flat liquid surface increases 
its surface area and takes a wavy form by inducing the diffusion of molecules from 
either bulk phase. Similarly a crystal surface is formed and reformed by processes of 
evaporation - by particle exchange with the surroundings. 

In fluctuating membranes such diffusion does not occur or is very rare. Com- 
monly membranes and bilayers are formed by self-assembling surfactant (amphiphilic) 
molecules and embedded in a liquid solvent. Although in principle the molecules forming 
a bilayer are soluble to some extent in the liquid solvent, in reality their concentration 
is exceedingly low. It is hardly possible for a bilayer or a membrane to absorb or release 
molecules in the process of locally changing its shape; consequently the shape fluc- 
tuations (also long-wavelength fluctuations, such as capillary waves) take place under 
the constraint of a constant particle number N s forming the sheet of surfactant. This 
constraint translates to a constraint of constant membrane area. 

Such considerations will make distinction between the true area and the "projected 
area" . The exact and formal results provided so far by the Gibbs statistical mechanics, 
always refer to the area defined in terms of the geometry of the macroscopic volume V 
enclosing the system under consideration. For example in d=3 the volume may be a 
parallelepiped and the area its x — y edge. A constant-volume increment 5F of the free 
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energy F provides the interfacial tension 7 through 

5F = 1 5A. (1.1) 

A is often called the "projected area", to emphasize that, because of the undulations 
of the interface, its "true area" is larger and, when projected onto the edge of the box, 
produces A. The "true interface area" is something the external world does not have 
control of; contrariwise, the projected area A is under the control of external forces just 
as the volume, temperature, or overall density is. Thus (1.1) defines 7 in terms of a 
measurable quantity A. Such reasoning clearly applies to an interface; we examine the 
other case of a membrane with a constant intrinsic area. 

These observations prompted us to consider a new set of SOS models with con- 
straints. In this paper we limit ourselves to one-dimensional string embedded in two 
dimensions, i.e. to a model which is exactly soluble[2-7] albeit without constaints. 

We solve below (for the first time) the microcanonical i.e. constant-energy version 
of the SOS model which we use to obtain the combinatorial factors; these are used to 
construct the new ensemble of constant interface length. Therefore we proceed via exact 
enumeration. 

In Section II the SOS model is defined and solved; this provides the combinatorial 
factor g(E; L, M) - the basic building block for subsequent calculations. 

These results/calculations have many ramifications, some of which are reported in 
Section II and in the Appendix because of their novelty and possible usefulness. 

In Section III the definitions of the total area are compared and then the relation 
between total area A to t and the projected area L is compared for the unconstrained 
calculation and for the new calculations under the constraint of a given constant Atot- 
Also, the derivative of the free energy with respect to either area (defining the tension 
7) is calculated and discussed. 

Finally Section IV is the summary and brief discussion. 
II. The Solid-On-Solid Model and its Solution. 

In the two-dimensional Solid-on-Solid (SOS) model[l-9] a one-dimensional interface 
is drawn along a strip of a square lattice as illustrated in Fig.l. The interface is reduced 
to a one-dimensional string. The standard L x M strip is made of L columns of height 



3 



M; the interface crosses each column once and only once. A microscopic configuration 
of the interface is then specified by a collection of heights 



{hj},j = l,L; l<hj<M. (2.1) 

Clusters and overhangs are absent. The "vertical" variable is h or z at the horizontal 
position x = aoj, 1 < j > < L and we take the lattice constant ao = 1. The "height" is 
limited to 1 < h < M. 

In the standard SOS model[2,5,6] the energy of a microscopic configuration is 

E tot = E + E (2.2) 

with 

E = e L; E = eA v . (2.3) 

In general e ^ e and both are positive. E is the energy of a perfectly flat interface. 

The total length of the interface is its one-dimensional area; it is a sum of its 
horizontal part, Ah and its vertical part, A v . 

Atot = A h + A v = L + A v (2.4) 

The first part is the trivial contribution equal to L. The vertical part is the sum of 
absolute values of the height jumps 

i i 

where 

A M+ i = \h i+1 -hi\ , (2.6) 

Thus the energy of the strip and of the interface is directly related to the vertical 
part of the interface A v . This is not essential. 

In more general SOS models where 

E = eJ2 \hi+i ~ H n e>0 (2.7) 
e.g. when n = 2 in the Gaussian model, the simple proportionality of E and A v , is lost. 
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Also, other definitions of the total area, can be given (see Section III) but are 
numerically close to (2.4). 

We define the following partition functions and their associated thermodynamic 
potentials. For the finite strip L columns long and M rows high, we have 

(A) the micro canonical constant-energy ensemble 

Z(L,M,E) = g(A v ;L,M) S = klogZ (2.8) 

where g is the number of microscopic configurations and S is the entropy, k is the 
Boltzmann constant, E = eA v . Given eo, one can use interchangeably Eq or L. Given 
e, one can use interchangeably E or A v . 

(B) The canonical (constant temperature) partition function is 

Z(T, L,M)= 9(A V ; L, M) x exp"^** (2.9) 

A v >0 

(3F = (-) log Z (2.10) 

where F is the free energy. This usually is written subtracting the energy of the flat 
interface e$L as the trivial term, defining 

Z(T, L,M)= 9( A ^ L ' M ) x exp"^ (2.11) 

a v >o 

so that F = —kT log Z is now the free energy of interface deformations. By introducing 
the quantity Q = exp[— f3e] G [0, 1] we rewrite (2.11) as 

Z(T,L,M)= g(A v )xQ A " (2.12) 

A v >0 

emphasizing that Z is the generating function for combinatorial factors g. Here, L being 
fixed, the total length of the interface A tot is directly given by A v and conversely. E, 
A v , and A to t fluctuate - appropriately for an interface. 

Generalizing the partition function for fixed L, we construct now 

(C) the new partition function for total interface area. That is, we consider an ensemble 
of strips with different lengths L (and common height M). Then 

Z*(T,A tot ;M) = ^^<7(A„;L,M)exp[-/3i? tot ] (2.13) 

L A v 
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with the implied constraint of a given total area, e.g. defined by (2.4) (as A tot = L+A v ). 
Total energy as given by (2.2)-(2.3) takes into account the energy cost of a horizontal 
step, e , besides the energy cost of a vertical step, e. It is convenient to define 



whereby exp(— (3E tot ) is written as Q Av R L . The partition function (2.13) becomes 
Z*(Q, R, A tot ; M). By allowing different energies for vertical and horizontal steps we 
can distinguish and treat as different, total lengths that are equal but are made of 
different proportions of vertical/horizontal steps. Both Q e [0, 1] and R G [0, 1]. Thus 
in (2.12) R = 1 has been tacitly assumed. The height M can be made infinite[7,8] 
for fixed L (Section lid). Under certain conditions, infinite M can also be introduced 
into (2.13). The averages are calculated as usual: < L >= —(dlogZ)f(dlogR) and 



To summarize, we consider the following geometries: 

(i) L fixed, M fixed and finite, 

(ii) L fixed, M infinite, 

(iii) L variable, M fixed and finite, 

(iv) L variable, M infinite. 

The new partition function for (iii) or (iv) with the constraint of given constant total 
area A to t = L + A V: is calculated in Section lid. 

In all cases we need the basic building blocks: the combinatorial factors g(A v ; L, M) 
which are computed in Section IIa,IIb, and g(A v ; L,M = oo) in Section lid. 

Ila. The Microcanonical Ensemble of Fixed Energy. 

To solve the SOS model in d = 2 for both M, L finite and given, one introduces the 
M x M transfer matrix as follows 



Q = exp[-/3e] R = exp[-/?e ], 



(2.14) 



<A V >=-(<91ogZ)/(<91ogQ). 



T(h,ti) =Q\ h ~ h '\ 



(2.15) 



with Q = exp[— e/kT]. The canonical partition function is expressed as [2] 




(2.16) 



h 
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If the transfer matrix is diagonalized, as a real symmetric positive definite matrix it 
admits a representation 

T hk = J2Mh)^Mk) (2-17) 

3 

where Xj is the j-th eigenvalue, j = 1, M and 4>j the corresponding eigenvector of the 
orthonormal and complete set. Therefore multiplication of T's just raises the power(s) 
of eigenvalues and 

Z = Af + A^ + --- . (2.18) 

The largest eigenvalue dominates in the limit L — > oo; then for sufficiently large 
M, F/(LkT) = -(l/L)logZ = -logAi = -log(l + Q) + log(l - Q) reproduces 
the known result, logtanh/3e/2 [2,7,10,12]. 

Now we want to select from the set of all configurations only those of given total 
length L + A v . If L is fixed this is reduced to a constraint of given A v . The latter is 

A v = A 12 + A 23 + --- + A Lil , (2.19) 

where 

A M+ i = \h i+1 - hi\ , (2.6) 

The last term in (2.19) is there if one imposes periodic boundary condition Jil+i = 
hi as we do. The constraint on A v can be expressed with the aid of a Kronecker delta 
5 Kr (j,n) equal to 1 for j = n and to otherwise. Thus 

Z(T, L, M; A V ) = J2-J2 A i2 + ...)* T hlth2 T T ■ . (2.20) 

h\ h L 

By using the integral representation 

{ ™ )e ikU-n) (221) 

and taking the integral in front of the sums, we find that to each T-matrix there corre- 
sponds a factor exp(i/cA) so defining the new matrix 

Th, h > = T Kh ,e lkA ^' , (2.22) 

we have the new expression for the constrained partition function 

(— )e"^Trr L . (2.23) 
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This is the expression we evaluate. The matrix r can be put into the form of the r.h.s. 
of (2.15); define Q as Q = Qexp(ik), so that 

r(h,h') =Q\ h - h '\ Q = e +lk ~^ . (2.24) 

The complex Q has a real part between -1 and 1 and an imaginary part between -1 and 
1. The matrix r depends parametrically on the wave- vector k, the integration variable, 
which is real. It reduces to the original T-matrix at k = 0, apparently continuously. 
By using the package "Eispack" [13], in particular its program CG for diagonalization 
of complex general matrices, we were able to compute Z and logZ from 

/ +7r Hk 
(^)e- lkA "[\i(k) L + \ 2 (k) L + --- + \ M (k) L } . (2.25) 

The eigenvalues depend on k. Evaluation was done for M = 5,6, 10, 20 and for L < 41; 
several values of Q were used. It could be done for M = 30. For odd A v the result 
is exactly nil and this is correct because of the periodic boundary condition at L. For 
A v = 0, Z = M for all Q. The general form of the result is 

Z(A V ) = g(A v )Q A " . (2.26) 

with g independent of Q. The Kronecker delta picks up the one term of the sum (2.12). 
By dividing (2.26) by Q Av we obtain g(A v ; L, M). All numerical results did conform to 
(2.26); g was obtained accurately for the very largest systems (g must be an integer) 
and was obtained exactly for all other sizes. Numerically, Q cannot be chosen too 
small or too close to unity. In actual calculations we have imposed periodic boundary 
condition h^+i = hi which not only simplified the calculation but also excluded exotic 
configurations. For all A v odd, g is nil. 

The combinatorial factor g(A v ) always follows the same pattern starting from 
g(0) — M, going through a maximum near A v ~ LM/3 and then falling down to for 
A v > (M — 1)L. Being the number of microscopic configurations, g very soon reaches 
horrendously large numbers and for that reason Fig. 2 shows an example of g(A v ) for 
very small M, L. The pattern seen in Fig. 2 is general; starting with g(0) = M and 
after a very rapid increase, g goes through a maximum and falls to a low value, often 
terminating with a value of 2. Afterwards for all values of A v the correct value is 
(7 = 0. This characteristic dependence becomes continuous for large L (and not too 
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small M ) in the following variables: Y = (1/L) log(g) and a = A v /Ln\. The quantity 
a = A v /Lni varies between and 1 , ri\ being equal to M — 1. 

Since the energy is not allowed to fluctuate, the logarithm of the partition function 
is the thermodynamic potential of the constant-energy micro canonical ensemble 

S(L,M,E) = klogg(A v ;L,M) E = eA v (2.27) 

An extension of the numerical results to larger systems, can be based on a simple 
scaling in reduced variables. We have already introduced the scaled, or reduced area 
(simply area per site) a = A V /L(M — 1) and we find that Y = (1/L) log(^) fall on a 
common curve when plotted against a, for different L at each constant M. 

Fig. 3 shows how this simple "scaling" works; the shapes are similar enough so that 
for each M different values of L fall on a common curve. Values of L < 10 and of M < 10 
are too small to be included. 

Eq.(2.27) identifies the function Y as entropy per column. Invoking the relation 

dS/dE = 1/T (2.28) 

defining the thermodynamic temperature T in the microcanonical ensemble, we may 
rewrite it as d\og(g)/dA v = f3e = — logQ with f3 = 1/kT. However, as a look at 
Fig. 2 and Fig. 3 shows, beyond the maximum in g that interpretation fails (or produces 
negative temperatures). The original observation is attributed to L.Landau[14]: in 
systems with energy bounded from above, T < is possible. Here is one more such 
example. 

The ordinary canonical partition function is calculated by diagonalizing the real 
T— matrix for given Q, i.e. after (2.16-18). But it is also equal to (2.12) . Since all 
terms are positive, selecting the maximum term leads to d(\ogg) / dA v — (3e = which 
can be rewritten to a form like (2.28). Now T is an independent variable, T > because 
Q G [0, 1], and E = E(T),S = S(T). 

lib. An Independent Algebraic Check. 

Turning to the starting point for the solution of the model, (2.15-2.16), we note that 
for small L the trace can be computed analytically. By multiplying out the T-matrices, 
we obtain a polynomial in Q: 

Z = TrT L = c + ciQ + c 2 Q 2 + ... . (2.29) 
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The basic observation is that the power of Q is just the number of energy units, identified 
in the previous subsection as the vertical part of the total length and denoted by A v . Z 
containing the total of all configurations, each term represents the split between different 
values of A v . That is , we can identify the coefficients Cj in (2.29) as the combinatorial 
factors denoted by g(A v ). Thus equation (2.29) is rewritten as 

Z = 9(A V )Q A " • (2.12) 

A„=0,l,2... 

For example, constructing the 4 by 4 T-matrix (whose first row is 1,Q, Q 2 ,Q 3 ) and 
multiplying it out to obtain T 6 , taking the trace, we obtain for L = 6, M = 4 Z = 
4 + 90x 2 + 510x 4 + 1266x 6 + 1116x 8 + 7Ux 10 + 310x 12 + 42x 14 + 12a; 16 + 2x 18 . Thus 
g(A v = 8) = 1116 etc. See Fig. 2. Such exact enumerations, with the computer-aided 
algebra, were carried out up to M = 20 for L up to 30, providing checks on numerical 
results obtained via the diagonalization described in previous Subsections. 

He. The Constraint of Total Area For Given M 

From the combinatorial factors in the micro canonical ensemble at given L, M we 
can construct now the partition function with the constraint of total length A tot . We 
take the ensemble of strips with different L's (and the same M) and implement the 
constraint of constant A tot in (2.13); for given M 

Z*(T, A tot - M) = J2/Z 6Kr( ~ L + A totMA v ; L, M)Q Av R L . (2.30) 

That is, the sum L + A v is held constant. Maximum value of L is obviously A tot and the 
mimimum value is 2. To avoid exotic cases, periodic boundary conditions are imposed 
tiL = hi. The Kronecker delta is replaced by these simple restrictions on the sums. 

Note that here we do not calculate the average total area A tot because the procedure 
is inverted: for given fixed total length A to t both L and A v are allowed to fluctuate and 
are calculated as averages: 

(L) = dlogZ* /dlogR (A v ) = d\ogZ*/d\ogQ (2.31) 

and (A v ) = A tot — (L) or conversely. The averages as f(Q,R,A tot ) were computed 
and are shown below for special cases such as u = Q = R and R = 1, Q e [0, 1]. In 
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the first case a string of fixed length is allowed to fluctuate in a neutral environment, 
where energy of a horizontal step is equal to that of a vertical step and therefore none 
is preferred; in the second case, only the energy of vertical steps is affecting the choice 
of configuration. 

Interestingly, the conbinatorial factors g(A to t) are increasing without bounds, in 
contrast to the constant L, M case from Section Ila-IIb. 

For an ensemble of narrow strips of M > 6 analytical calculations were done as 
follows: first traces of T L were calculated and stored (for L up to 40 and nore) and then 
sums were formed for all even values of Atot £ [10,38]. From the resulting expression 
j3F, (L), (A v ) were calculated analytically and used, in part to check the results of nu- 
merical computations. A selection of these results is shown in Section III. The results 
for finite M are also compared with the case M = oo described below. 

lid. The Columnar Strip (M=oo) and its Solution. 

From the solution [7] of the eigenvalue problem of T for given L and columns of 
infinite height [7], M = oo, the eigenvalues are known explicitely 

X(u) = (1 - Q 2 )/(l + Q 2 -2Q cos(i/)) (2.32) 

where the index v e [— it, +7t] is now continuous. Because of the translational invariance 
in the "vertical" M direction, actually exploited in solving the transfer matrix[7], we 
must impose a constant hi, e.g. hi = and another fixed value for h^ - which we choose 
to be equal to hi. The canonical partition function (2.11) or (2.12) with M infinite then 
exists. There are some limitations in more general cases (R ^ 1 and variable L fixed 
A-tot) - because now the excursions of the interface are limited only by the cost of a 
vertical step. The flat configuration hi = is the equilibrium configuration for Q — > 0. 

We need the combinatorial factors g(A v ; L, M = oo). Following Section Ha, in order 
to impose the the constraint of fixed A v we replace Q = exp(— (3e) with Q = Qe tk . The 
solution[7] goes through as explicit calculation shows. The partition function becomes 
a double integral, replacing (2.25) by 

Q)e~ ikA " J {^)m L ■ (2-33) 
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Numerical evaluation is now simpler as the diagonalization step is not needed - we have 
the explicit expression (2.32) for the eigenvalues. For A v and L integers, we obtain 
again Z in the form (2.26) i.e. Z = g x Q Av , where again g(A v ) comes out an integer. 
Again logg has the interpretation of entropy S/k. 

In an analytic calculation, all g(A v ;L) result as polynomials in A v . There is a 
striking simplification which occurs for M = oo; g(A v ; L) are polynomials with highest 
power equal to L — 2. For example g(x; L = 4) = 2 + (2/5)x 2 (for all x > 0), g(x; L = 
6) = 2 + (35/8)x 2 + (21/32)x 4 etc. (x = A v ). Knowing these polynomials we can make 
computations with any value of A v . A few examples are given in the Appendix. 

The exact results for g(A v ; L, M = oo) are used now to construct an ensemble of 
columnar strips of infinite heights M and variable L. At given Atot, maximum value 
of L is Atot, and the mimimum value is 2. Minimum value of A v is zero; there is no 
upper bound other than A tot — L. To avoid exotic cases and to ensure the existence of 
Z, again Jil = hi = for each L. We show here, in Section III, explicit calculations for 
these two cases: Q = R — > u and R = 1, Q G [0, 1]. 

III. Numerical Results and the Comparison of Areas and Free Energies. 

In this section we show a restricted selection of results for the two most important 
quantities: the vertical area A v related to the energy and the total area of the interface. 
We also consider the derivative of the free energy F with respect to area, all other 
variables being kept constant; this is interpreted as the interfacial tension 7. Hence 

07 = (d(pF)/dA) T (3.1) 

The area A can be the projected area L, the total area imposed in the ensemble of all 
L's, or the average total area defined as (A tot ) = L + (A v ). All results are given for 
R = 1,Q e [0,1] (cf.(2.13)-(2.15)). The other interesting case, R = Q E [0,1], turns 
out to be a degenerate temperature-independent one, in which the interface wiggles 
randomly in an energetically neutral environment. These results are also reached in the 
limit Q — > 1 _ at constant R = 1. 

First, however, we compare 3 different definitions of the area, or, rather, two def- 
initions, A v and A sq , and one approximation, A' . Our definition of the total area as 
A to t = L + A v corresponds to the length of the thick line in Fig.l; L is the sum of 
"horizontal" steps, and A v is the sum of "vertical" steps, A^+i (cf. (2.5-6). The usual 
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definition would draw a line through hi and fti+i; this approach produces 



A *i = + y/ 1 + a m+i]- 



(3.2) 



where L is subtracted. Expansion of the square root gives the first term, which is the 
common approximation, 



Definitions A v and A sq produce always very close results, A v always slightly larger. 
Fig. 4 shows all three areas calculated in the canonical ensemble, plotted against Q. 
A' is a very bad approximation, grossly exaggerating A sq except when both are very 
small; this happens here at small Q. We have found earlier in another context [15] the 
same pattern, with the commonly used approximation for the area, like (3.3), grossly 
overestimating the true area, except in the limit AA — > 0. 

Figure 4 also includes the microcanonical results i. e. the curve of A v as independent 
variable against Q = exp(— {j3e)) with the inverse temperature computed acccording to 
(2.28) with (2.27); even for the small value of M the differences are negligible and 
invisible on the scale of the plot. This is a useful check; the agreement between the 
canonical and microcanonical ensembles can only become better for larger L and M 
and it does. In conclusion, we use A v + L as a measure of area, include canonical 
averages, and do not quote microcanonical averages. 

Fig. 5 illustrates the results for the new ensemble of variable L and constant A tot . 
Three values A tot = 10, 20, 30 were chosen. Since A sq and A v are so close, only (A v ) is 
plotted. The case M = oo (section lid) is shown with lines; at the infinite temperature 
(A v ) reaches 4.2, 9.2, and 14.2 for A to t = 10, 20, 30, respectively. Interestingly, squeezing 
the membrane with a small M = 6 does not change (A v ) much; these averages are shown 
with stars, crosses, and plus signs, and reach 4.0, 8.6, and 13.0. Unconstrained averages 
calculated in the canonical ensemble produce plots of another shape and much larger 
values at high temperatures; (A v ) for M = oo diverges. For fixed A to t it does not. 

We choose the inverse temperature as the ordinate for the plot because of the 
quasi-linearity at high temperatures of the data with fixed Atot- As can be seen clearly, 
fixing A tot limits the vertical excursions of the interface considerably and - somewhat 
suprisingly - in a similar fashion for small M = 6 and for M = oo. At low temperatures 
i.e. at vanishing Q and large (3e all calculations produce (A v ) — > 0. 




(3.3) 
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At Q = 1 i.e. (3e = 0, (L) is at 55-57 percent of A totl practically in all calculations 
we have done so far for the new ensemble. The limit corresponds to the fluctuation 
of the interface in the neutral environment where horizontal or vertical steps carry the 
same cost in energy, now vanishing in the limit Q — > 1, j3e — > 0. 

The divergence of the canonical average for M = oo (full line in Fig. 5) can be 
understood as follows. If M = oo and L is fixed whereas A v and A tot are free to 
fluctuate and take an average value, it is only the energy cost expressed in the value 
of Q which limits the excursions of the interface and therefore the number of realized 
configurations thus the partition function. When Q ceases to operate, the partition 
function and the averages, diverge. This does not happen when a fixed A to t is imposed. 

All partition functions are converted to free energies and all free energies depend 
on an area as the independent extensive variable. The definition (3.1) of the interfacial 
tension 7, is now applied to (3F(A tot ;T). In Fig. 6 the derivative (3.1) w.r.to A tot is 
plotted against Q. For M = 00 or M = 6 in the ensemble of all L, the data points show 
mild near linear variation with Q. 

For comparison, the canonical free energy per column [3F(Q; L, M) / L from the 
largest eigenvalue, is — log[(l + Q)/(l — Q)] and is plotted as the dot-dash line. All 
eigenvalues for M = 6 produce the near-linear full line. 

Whether the derivative (3.1) applied to F(A tot , T) can be interpreted as an "interfa- 
cial tension" , is perhaps debatable. Still the derivative can be calculated and compared 
with other derivatives of the free energy with respect to other areas. 

IV. Summary and Discussion. 

As recalled in the introduction membranes keep their area constant, notwithstand- 
ing shape fluctuations. Planar interfaces are also twodimensional sheets embedded in 
three dimensions, but are open to particle exchange with the surroundings and their 
proper or intrinsic area varies. Hence the theory of fluctuating interfaces must be suit- 
ably modified to apply to membranes. We examine this modification by choosing an 
existing simple and soluble model, namely the SOS model in two dimensions, for which 
exact enumerations were practical. 

By starting from the available results, we introduced the constraint of constant 
interface area; first simply into the L x M strip and then constructing a new ensemble 
of all projected areas. 
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The partition functions, free energies, and the ensemble averages all were built on 
the microcanonical combinatorial factors g(A v ; L, M). Their calculation was essential. 
These were obtained by several means described in Section II. Of the several calcula- 
tional schemes, the ensemble of all lengths L at infinite height M under the constraint 
of a given total area A to t produced (A v ) and (L). The latter has an interpretation of the 
average projected area. Practically in all calculations, the ratio (L)/A tot approached 55 
percent in the limit of high temperatures (Q — > 1, f3e — > 0). 

We also calculated the derivative of the free energy F(A tot ,T) with respect to 
A to t- This may be considered an analogue of the interfacial tension 7, now defined with 
respect to A to t- When compared with other calculations of similar derivatives, the new 
derivative was the smallest in the absolute value and did not diverge in the limit of 
infinite temperature. This was illustrated by Fig. 6. 

The mathematical procedures involved some particular features; their description 
may be useful. First, we used a diagonalization of a complex matrix r which was 
symmetric but not unitary. This was followed by numerical integration. The final result 
had to be an integer. The correcteness of the diagonalization procedure was verified 
by agreement with the exact algebraic calculations. The exact enumerations (for the 
case of infinite M) were also spot-checked by the generating functions for "partially 
directed" random walks, obtained by slightly modifying the generating functions given 
by Privman and Svrakic[18]; but introduction of constraints was, as often is the case, 
not any more practical. 

Extension to two-dimensional membranes embedded in d = 3 is not impossible 
though difficult. A good candidate[17] is the Gaussian SOS model on a square lattice 
with each site endowed with a continuous height variable hf, it is readily soluble without 
constraints and for various boundaries. The prediction is that the constraint will damp 
the fluctuations in a way similar to the damping seen here. 

Finally, it is probably all too obvious that our aim was not to develop further the 
Solid-on-Solid model but rather to use it to examine the differences between membranes 
and interfaces. Thus we have not pushed the calculations in order to go to highest 
possible orders of calculation . The trends are already visible. 

Acknowledgements. 

The author has greatly benefitted in the past from discussions with J.S.Rowlinson 



15 



FRS (Oxford), R.Evans FRS (Bristol), and B. Widom (Cornell). 



Appendix. 

Combinatorial Factors g for Infinite M. 

It is a rare and desirable event if combinatorial factors can be generated by polyno- 
mials. We list the first few combinatorial factors g(x) for L =2, 4, 6, 8, 10. These were 
calculated and are available for L up to 30. These polynomials are valid for any value 
of x = A v > 0. For x = A v = 0, g = 1 for any L. Some exact algebraic calculations 
were done in part with the aid of an old version of " Mathematica" [16]. See Section lie 
for notation and further details. 

L= 2 g = 2 

L= 4 g = 2 + (5/2)x 2 

L= 6 = 2 + (35/8)x 2 + (21/32)x 4 

L= 8 = 2 + (707/120)x 2 + (77/48)x 4 + (143/1920)o: 6 

L=10 g = 2 + (14465x 2 )/2016 + (4147x 4 )/1536 + (715x 6 )/3072 + (2431x 8 )/516096 ... 
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Figure Captions 
Figure 1 

The solid-on-solid (SOS) interface in d = 2. The strip is M sites wide in the 
"vertical" direction and L sites long in the "horizontal" x direction. The interface is of 
length A = L + A v . 

Figure 2 

An example of the number g of configurations of the interface plotted against the 
interfacial area A v at given size of the strip with L=6,M=4. The numbers (different 
from zero) are: 4, 90, 510, 1266, 1116, 744, 310, 42, 12, 2, for areas A v 0,2,4,6,8,. ..,18, 
respectively. 

Figure 3 

The "scaling" curve of the microcanonical combinatorial factor g = g(A v ) at fixed 
L,M. The function Y = (l/L)log(#) is plotted against a = A V /L(M - 1) for M = 
6, 10, 20, and L = 10, 20, 30. Smaller values of L deviate from the common curve. The 
maxima extrapolated to 1/L — > with polynomials in 1/L were 1.77702, 2.28587, and 
2.97979 for M = 6, 10, 20, respectively. Extrapolation to 1/M — > was inconclusive. 

Figure 4 

Three areas: A v /L, A sq /L, and A' / L vs. Q calculated in the canonical ensemble 
(see text) are compared. Full line terminating at 5.8198 - A sq /L; bold-dashed line 
hitting 33.25 - A'/L; fine-dashed line hitting 6.6500 - A v . Plus signs - microcanonical 
A v reaching 6.564 at Q = 1. L = 30, M = 20. A v /L was independent of L and weakly 
depending on M provided M > 10. 

Figure 5 

Average (A v ) vs. f3e = e/kT = — logQ for fixed Atot = 10, 20, 30 in the ensemble 
of all L. Fixed M = 6 shown with points (stars, crosses, plus signs, respectively); M = 
oo shown with lines (fine-dashed, dashed, full, respectively). Unconstrained canonical 
averages (A v ) for constant L, M are (M = 6, L = 10, 20, 30, - dotted lines; ditto M = 
oo, L = 10, - full line) . For /3e = 0, data shown with points hit 4.00, 8.564, 12.980, 
respectively; data for M = oo shown with lines, hit 4.205, 9.209, and 14.213, respectively. 
The dotted lines (canonical (A v ) for M = 6) hit 19.444, 38.888, 58.333 for L=10,20,30, 
respectively. R = 1, Q G [0, 1] in all cases. 
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Figure 6 

Free energy derivatives w.r.to the total area (d/3F/dA tot )Q, plotted against Q. The 
new ensemble of all L and M = 6 - small plus signs, ditto M = oo -large plus signs. For 
comparison the derivatives w.r.to the projected area L, in canonical ensemble M = 6 - 
full line, M, L — > oo - dash-dot line, are plotted. 

the end of text. Figures are as 6 *.eps files 
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